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Abstract 

A disclination in a hexatic membrane favors the development of Gaussian cur- 
vature localized near its core. The resulting global structure of the membrane 
has mean curvature, which is disfavored by curvature energy. Thus a mem- 
brane with an isolated disclination undergoes a buckling transition from a flat 
to a buckled state as the ratio k/ Ka of the bending rigidity k to the hexatic 
rigidity Ka is decreased. In this paper we calculate the buckling transition 
and the energy of both a positive and a negative disclination. A negative 
disclination has a larger energy and a smaller critical value of h/Ka at buck- 
ling than does a positive disclination. We use our results to obtain a crude 
estimate of the Kosterlitz-Thouless transition temperature in a membrane. 
This estimate is higher than the transition temperature recently obtained by 

the authors in a renormalization calculation. 
PACS numbers: 05.70.Jk, 68.10.-m, 87.22.Bt 
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I. INTRODUCTION 



The hexatic phase [0] is characterized by 6-fold orientational but not translational order. 
Both three-dimensional hexatic phases with true long-range order and two-dimensional 
phases with power-law order have been observed. In flat two-dimensional films (under 
tension), the transition to the isotropic phase occurs via a Kosterlitz-Thouless (KT) discli- 
nation unbinding transition |^. Free membranes with zero surface tension can also exhibit 
a hexatic phase and a KT transition to a fluid phase. Both the hexatic phase and the tran- 
sition to the isotropic phase are, however, more complicated than they are in a fiat film 
because of thermally induced shape fluctuations. Recent renormalization group calculations 
show that shape fluctuations shift the bare hexatic rigidity Ka- As a consequence, 
an increase in the amplituded of shape fluctuations produced by decreasing the membrane 
bending rigidity /t will induce a transition from the hexatic phase to the isotropic fluid phase. 

In flat membranes, there is a symmetry between positive (5-fold) and negative (7-fold) 
disclinations, and they both have the same energy. In free deformable membranes, this sym- 
metry is broken ||^. If the ratio h/Ka of the bending rigidity k to the hexatic rigidity Ka 
is sufficiently small, a membrane with a single dislcination can lower its energy by buckling, 
thereby creating a nonvanishing Gaussian curvature, which screens the disclination charge. 
The buckled states of positive and negative disclinations have different height profiles, dif- 
ferent energies, and different critical values of n/ Ka at which buckling occurs. In this paper, 
we will use a variational procedure to calculate the energies of isolated positive and negative 
disclinations on free membranes. Our variational form for a positive disclination is essen- 
tially exact. Our form for negative disclinations is exact only for k/Ka near the critical value 
of njKA for buckling. Corrections near buckling are, however, very small, and we argue that 
our form is a very good approximation until k/Ka becomes very small. Our results are in 



agreement with recent calculations by Deem and Nelson |T0|. The latter authors, however, 
in addition to calculating the energy of a negative disclination variationally also carry out a 
numerical minimization of the full non-linear energy to obtain a lower value for this energy 
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at small values of h/Ka- In what follows, we will present our calculations of positive and 
negative disclinations, respectively in Sees. II and III. In Sec. IV, we will discuss our results 
and an estimate of the KT transition line produced by them. 

The continuum Hamiltonian for hexatic membranes was derived in Ref. [0]. If we 
parametrize membranes positions in terms of a two-dimensional parameter u = 



as R(m) , then 



where 



n = n. + nc, (i.i) 



n. = l^^J d\^H\ (1.2) 



is the curvature energy and 



nc = -KA J d'u^{^-S)-^{^-S), (1.3) 



is the Coulomb energy. Alternative form for the hexatic energy can be found in Ref. |]TT 
In the above, g = det{gab) is the determinant of the metric tensor gab = ■ i9bR , 
Af = ^nj^ili^iu ~ '^i)! \f9 is the disclination density with = ±1/6, and H and S are, 
respectively, the mean and Gaussian curvatures of the membrane. We have ignored a scalar 



field contribution of 7i, which gives rise to Liouville measure factors [jl2|,|T3[, which are 
irrelevant to the current discussion. On a rigid flat membrane, this Hamiltonian reduces 
to the Coulomb gas form of the XY-model. In the absence of disclinations, hexatic order 
induces long-range Coulombic interaction between Gaussian curvatures on the membrane 
0. The Coulomb energy Tie depends only on the difference M — S. Thus, the development 
of Gaussian curvature on a free membrane that approximates the disclination dentsiy M 
can reduce the Coulomb energy. Gaussian curvature usually leads to mean curvature, and 
the lowest energy state of a free membrane with a disclination will be determined by the 
competition between the Coulomb energy. Tie*, which prefers S = JV, and the curvarure 
energy, Ti^, which prefers zero curvature. If there is a single disclination at the origin, one 



can expect that Gaussian curvature will be localized near the origin. Gaussian curvature 
localized in a small region will give rise to a buckled state with mean curvature but zero 
Gaussian curvature away from the origin. 



II. POSITIVE DISCLINATIONS 

A minimum strength positive disclination has "charge" g = 1/6. To reduce the Coulomb 
energy associated with this charge, the membrane can distort into the shape of a spherical 
section, with nonvanishing Gaussian curvature, localized to the core of the disclination. 
Outside the core region, the membrane will seek a shape with zero Gaussian curvature. A 
cone with slope m has zero Gaussian curvature and can be connected smoothly to a spherical 
section (Fig. |I]). Thus, in the Monge gauge, we parametrize the membrane shape outside 
the core as -u = (r, 0) and R('u) = (r cos0,rsin0, /i('u)) with h{u) = mr. Thus, outside the 
core, the components of metric tensor and its determinant are 

gre = 9er = 0, g = r^{l + m^). 

The mean curvature H and the Gaussian curvature S are 

g= , -, ^ = 0. (2.2) 

Thus the bending energy of the cone with radius R and the core size a becomes 

1 27rrdr , ( m \ 

= -K V 1 + ^ 

Ja 



2 Ja r2 \Vl + my 

w? ^ R , , 

= vr«:^==ln-. (2.3) 

The Gaussian curvature vanishes outside the core region. In the limit of the infinitesimal 
size of the core region, S can be described by the point curvature charge s+: 

S{u) = 2ns+5{u - u+)/^/g. (2.4) 



Since there is no Gaussian curvature in the cone, we can choose any curve in the cone to 
calculate s+ using the Gauss-Bonnet theorem: 

/ Sda+ [ Kdl = 2n, (2.5) 
Jm Jc 



where kg is the geodesic curvature of the boundary curve C of the surface M |T^. We use 
the boundary curve of the cone: 

C = (i? cos sin 0, mi?). (2.6) 



The geodesic curvature of the boundary curve of a cone of slope m is {Ryl + m?) ^, and 
the Gauss-Bonnet theorem becomes 



/ Sda+ , = 271. (2.7) 

Jm Jo VI + m2 

Using Eq. (p^.4|), we obtain 



V 1 + ^ 



Hence, the Coulomb energy for a positive disclination is 

He = \ka j d^u^{M - S)^^{N - S) 



7rir^(^-s+y27rGe(0), (2.9) 



where s+ = 1 — l/vT+m^ and Gc{u) is the Green's function for the Laplacian, 

V' = ^d^g''^^df,. (2.10) 
V 9 

On a cone, 

V^Gc{u-u') = -S{u-u')/^. (2.11) 

To determine Gc, we assume that it has the form —Aln{r /ro) where tq is a length. Then 

J d\y^V'G,{u) = - j d\^6{u)/^ = -1 

= J dsag^'VadbGc =-j dsr-g^'A/r 

' dSr{geel^){A/r), (2.12) 



where dsa = Sarrd6 is the "surface" element of a circle enclosing the orgin. Then from Eq. 



Q, gee/^ = r/VT+^, A = -VT+^/{27c), and 

GJu) - 



27r 



In In — 

a a 



(2.13) 



where we chose tq to be equal to the disclination core radius a and we added the constant 
term — ln(i?/a), where R is the radius of the cone, to produce the required divergence of 



Gc{u) at small r. The Coulombself-energy is given in terms of 27rG'c(0) = vT+m?ln(i?/a) 
where R is the radius of the membrane and a is the core size. 

The energy of a positive disclination with q = 1/6 on the cone with the slope m becomes 

.2 /I , X2 



Ejm) 



m 



7TK 



+ ttKa - - 1 + 



1 



6 



1 / K 11 

36 ^ ~ 72 



VTT 

83 



VT 



1 K 



I m 



In 



a 



R 



(2.14) 



.288 2KaJ 

This energy is shown in Fig. ^(a) for various values of k/Ka- For k./Ka > 11/72, Ej^{m) has 
a minimum at = and the membrane remains fiat. For k/Ka < 11/72, however, Ej^{m) 
has a minimum at m? = (11/36 — 2k,/ Ka)/ {25/36 + k/Ka) and the membrane buckles out 
to form a cone with the slope 



n 

36 



2k 
Ka 



(2.15) 



-\ 25 I K ■ 
\ 36 ^ Ka 

Thus the buckling transition occurs at k/Ka = 11/72 for positive disclination with q = 1/6. 
This result for m with Eq. p.l4| can be found in Ref. [jl5| and [|l^. The energy of a positive 
disclination on a membrane of radius R with the short- distance cutoff a is 



36k 
25^:^ 



^ttKaIu^, 

36 ^ a ' 



1 ) In^, 

a ' Ka 72 

K ^ 11 

Ka 72- 



(2.16) 



When Ka — > oo, m ±yll/25 and s+ 1/6. Thus, in this limit, the disclination 
charge is totally screened by the Gaussian curvature, there is no Coulomb energy, and the 
disclination energy Ej^{Ka = oo) = (ll/30)7rK ln(i?/a) comes entirely from curvature of the 
cone. This energy of a positive disclination is identical to the result obtained by Guitter 
and Kardar using the conformal gauge. 



Our simple height profile for a positive disclination does not break azimuthal symmetry, 
and there is no particular reason for this symmetry to be broken. Thus, we believe that 
h{u) = mr provides a complete discription of the buckled state and our description of the 
positive disclinations is exact. In particular, no symmetry breaking terms such as rriir cos 20 
or cos 40 are needed in the expansion of h{u). Order parameters such as mi and m2 are 
certainly not forced by the development of non-zero m because symmetry does not permit 
terms linear in of the form m^mp to appear in the expansion of E^{m). 

A KT melting temperatue T+ for positive disclinations can be introduced in the usual 
way by setting the free energy of a single disclination equal to zero; — = 0, where 
S = ln(_R/a)^ is the entropy. Thus = E+/2\n{R/a). This produces the phase diagram 
obtained in Ref. and shown as the solid curve (a) in Fig. ^. The thin line indicates 
the buckling transitions at k/Ka = 11/72, and the solid curve the disclinations unbinding 
transtion obtained from T+. 

III. NEGATIVE DISCLINATIONS 

We can similarly calculate the energy of negative disclinations with g = — 1/6. Since the 
corresponding core region should have a negative curvature charge to cancel the topological 
charge, we expect the core region has a saddle shape. The simplest saddle shape (Fig. ^) is 

h{u) = mr cos 2(p . (3.1) 

We will take this as a variational function and seek the minimum energy solution for a 
negative disclination with respect to variations in the parameter m. We thus obtain an 
upper bound to the energy of a negative disclination. Inclusion of additional terms in 
h{u) proportional to cos2n0 for n an integer will lead to lower energies. Indeed, recent 
numerical calculations by Deem and Nelson [|T0[ yield a lower energy than we obtain when 
Ka/k, ^ 1. We will argue, however, that this simple variational form is essentially exact 
near the buckling transition. 
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The components of the metric tensor and its determinant associated with h{u) — 
mr cos 20 are 



Qrr — i + vr? cos^ 20, geQ — r^{\ + rvr? sin^ 20), 

9rd = 9er — —2rm^ COS 20 sin 20, g = r^(l + m^(l + 3sin^ 20)). 

The mean curvature H and the Gaussian curvature S are 

1 + cos^ 20 



(3.2) 



_ 3m 

" -— cos 20^^ ^ ^2(1 + 3sin2 20))3/2 



The bending energy of the saddle with slope m is 



, ,5 = 0. 



(3.3) 



1 r^rdr r^^ Qm^ cos^ 2(I){1 + cos^ 2(l)y 
^2" la T^io '^'^ (l + m^d 



(l + m2(l + 3sin2 20))V2 



9 2 
2 



cos^ 20(1 + cos^ 20)^ 



a 



(3.4) 



/o ^(l + m2(l + 3sin2 20))V2_ 

Again, in the limit of the infinitesimal size of the core region, S can be described by the 
point curvature charge s_: 



S{u) — 2tts_S{u — uS)l 



(3.5) 



The integrated geodesic curvature along the boundary C — (i? cos 0, sin 0, mi? cos 20) is 

(1 + 4m2)rf0 



\ kgdl — -r- 

Jc ^ Jo 1 



/C - Jo (1 + m2cOs2 20)1/2(1 + 4^2 sin2 20)3/2- 

Thus the Gauss-Bonnet theorem gives 

1 /■2T (l + 4m2)d0 



(3.6) 



s_ = 1 



2^ Jo (T 



(1 + m2 cos2 20)1/2(1 + 4^2 sin^ 20)^/2 ' 
and the Coulomb self-energy of negative disclination on the saddle becomes 



(3.7) 



He = Ik A I d\^{N -S)^{N-S) 



7^Ka[- + s_] 27rG,(0), 



(3.8) 
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where s_ is given by Eq. ( |3.7| ) and Gs{u) is the Green's function for the Laplacian on the 
saddle. We can determine Gs{u) following exactly the same procedure we used to determine 
Gc{u) for a positive disclination. We find 



G,{u) = -A[ln{r/a) - \n{R/a)] 



(3.9) 



where 



A 



1 + 4m^ sin 



1 

2^ 



"^[1 + m2(l + 3sin2 20)] V2^ 

3 2 67 4 
1 m H m + ■ ■ 

4 64 



(3.10) 



Equations ( p.4| ), (|37^ ), (|3.9| ), and ( |3.10|) completely determine the energy of a negative 
discliantion as a function of m within the approximation h{u) = mr cos 20. We can locate 
the buckling instability and determine m just above it by expanding E^{m) in powers of m 
up to order m^. The result is 



1 ^ 9 ^ K 13 2 ^ (^2743 207 k 



m 



In 



R 



(3.11) 



.36 2 \Ka 216/ ' V2306 16 Ka. 

This energy is shown in Fig. |^(b) for various values of h/Ka- For k/Ka > 13/216, E^{m) 
has a minimum at = 0; for k/ Ka < 13/216, E^{m) has a minimum at m? = ((13/216) — 
/t/ir^)/((2743/10368) - (23/4) ■ k/Ka). The buckling transition occurs at 



K _ 13 
^ ~ 216' 



(3.12) 



and the slope for k/Ka < 13/216 is 



(13/216) - k/Ka 



(3.13) 



(2743/10368) - (23/4) ■ k/Ka 

The energy of negative disclination around k/Ka = 13/216 in a membrane of radius R 
with the short- distance cutoff a is 



E^ 



^ (l 9((13/48)-(9/2)(k/Ka))" \ 1 R _K_^^o/r,,a 

36-^^^ V ((2743/2304)-(207/16){K/_ft:A))y a' Ka ' 

5k^Alnf, ^> 13/216. 



(3.14) 



As in the case of positive disclinations, the Gaussian curvature will adjust to exactly cancel 
the topological charge when Ka = oo, leaving only curvature energy. Setting s_(m) in Eq. 
( p.TD equal to —1/6, we obtain 

m^{KA = 00) = 0.350417 (3.15) 

and 

R 

= 00) = 2.75883k In-. (3.16) 

E^{m) can be minimized numerically for < k/Ka < 13/216. The results are displayed as 
a transition temperature in Fig. |^ (See below). 

As discussed in the introduction, the height profile of a negative disclination breaks 
azimuthal symmetry, and we expect h{u) to have a Fourier series expansion of the form 
h{u) = r J2n^n cos 2n6. Our approximation keeps only the first term in this series. Near 
the transition, higher order terms can be calculated by expanding E_ in a powers series in all 
of the m„'s. We have already calculated the contribution from the dominant term mi = m. 
One might expect that the next most important term would be m2- This parameter is not, 
however, forced to develop a nonzero value when m is nonzero because E_ is invariant under 
h —h, i.e., under m„ — — m„ for every n. Thus, there is a contribution to of the form 
a2m| but no term proportional to m^m2, which would force a nonzero 1712. Thus 1712 will 
remain zero until the coefficient a2 changes sign. The absence of an m^m2 term means that 
our expressions for E^ [Eq. ( |3.14|) ] and m [Eq. (|3.13| )] are exact to order [(/t/i^^) — (13/216)]^ 



because there is no correction to the term arising from couplings to m2. The third order 
term is proportional to because inversion symmetry permit a term proportional to 
m?m^. More generally, there are couplings of the fomr m'^'P^^m2p+i for p = 1,2, ... . Thus 
the height profile can be expanded as h{r, 6) = r I]p=o "^2p+i cos[2(2p + 1)6']. Our results for 



E^ and m agree with those obtained analytically and numerically by Deem and Nelson [jlO 
Their numerical result for n/ Ka = is lower than our indicating that the order parameters 
m2p+i for p > 1 are important in this limit. 
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A KT transition temperature for negative discliantion can be introduced just as for 
positive disclinations: T_ = E_/2\yi{RI a). This curve is shown as the sohd curve (b) in 
Fig. |. 

IV. DISCUSSION 

Clearly both positive and negative disclinations will be thermally excited, and neither 
T+ nor T_ is a good estimate of the actual melting temperature, Tm- A better estimate is 
that Tm is simply the average (T+ + T„)/2. This yields the dashed curve in Fig. |^. This 
estimate describes qualitatively features that are in agreement with simple physical reason- 
ing: for large k, there should be a disclination-mediated melting to the disordered phase 
as temperature is increased, and at fixed Ka, there should be a transition to the crumpled 
phase as k is decreased. In Refs. |0] and @], we calculated the melting temperature using 
the renormalization group (RG) recursion relations for the KT transition on a fluctuating 
membrane subject to the constraint of charge neutrality. The result shown as the solid curve 
in Fig. 1^ is below the estimate (T+ +T_)/2 in the entire region below the dotted curve where 
we believe that our RG calculations are valid. This is entirely reasonable. The estimate of 
Tm obtained by equating the energy of disclinations to temperature times their entropy 
completely ignores the entropy associated with height fluctuations, which should lead to a 
depression of Tm- Our RG calculations include height fluctuations, whose major effect is to 
decrease the effective long- wavelength dielectric constant. 
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and negative disclination energies and to Michael Deem and David Nelson for sending us 
their independent work on this subject prior to publication. This work was supported in part 
by the National Science Foundation under grant No. DMR94-23114 and in part by the Penn 
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FIGURES 

FIG. 1. Membrane buckled into a cone with h{u) = mr in the Monge gauge where m is the 
slope. 




(a) (b) 
FIG. 2. (a) Energy of positive disclinations as a function of for different values of p = k/Ka'- 

(1) p/pt = 4/3, (2) p/pt = 1, (3) p/pt = 3/4, and (4) p/p+ = 1/2, where pt = 11/72. (b) Energy 

of negative disclinations as a function of m? for different values oi p = k/Ka'- (1) p/ Pc = 4/3, (2) 

p/p- = 1, (3) p/p- = 3/4, and (4) p/p- = 1/2, where p- = 13/216. 
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0.01 0.02 0.03 0.04 0.05 
FIG. 3. Estimated phase diagrams in the (T/k,T/Ka) plane showing the Kosterlitz-Thouless 

transition line obtained by balancing energy and entropy of (a) a single positive disclination 
(T+ = £'+/21n(i?/a)) and (b) a single negative disclination (T_ = E^/2ln{R/a)). (a) is identical 
to the estimate obtained in Ref. ||6| . The straight line through the origin in both cases is the buck- 
ling transition line. The energy of a negative disclination is generally higher than that of a positive 
disclination, and r_ > T+. T+(Ka = oo) = (ll/60)7r/t ~ 0.575959k and T^{Ka = oo) ~ 1.37941k. 



14 



FIG. 4. Membrane buckled into a saddle with h{u) = mr cos2(f> in the Monge gauge where m 
is the slope. 



r/K 




0.01 0.02 0.03 0.04 0.05 

FIG. 5. Phase diagram in (T/ k,T/ Ka) plane. The full line is the KT transition line obtained 

from the RG calculation of Refs. and |^. The dashed line is the estimate of the KT transition 
line from the average (T+ + T_ ) /2 of the positive and negative melting temperatures T-f and T-_ . 
The approximations used in Refs. and § apply below the dotted line. The RG transition line 
(full line) obtained in Refs. and Q lies below the simple estimate (dashed line) in the region 
(below the dotted line) where the RG calculation applies. This is expected since the RG calculation 
includes contributions to the free energy arising from thermal membrane fluctuations which the 
simple estimate (T+ + T_)/2 does not include. 
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